Breakdown of the few-level approximation in collective systems 
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The validity of the few-level approximation in dipole-dipole interacting collective systems is dis- 
cussed. As example system, we study the archetype case of two dipole-dipole interacting atoms, 
each modelled by two complete sets of angular momentum multiplets. We establish the breakdown 
of the few-level approximation by first proving the intuitive result that the dipole-dipole induced 
energy shifts between collective two-atom states depend on the length of the vector connecting the 
atoms, but not on its orientation, if complete and degenerate multiplets are considered. A care- 
ful analysis of our findings reveals that the simplification of the atomic level scheme by artificially 
omitting Zeeman sublevels in a few-level approximation generally leads to incorrect predictions. 
We find that this breakdown can be traced back to the dipole-dipole coupling of transitions with 
orthogonal dipole moments. Our interpretation enables us to identify special geometries in which 
partial few-level approximations to two- or three-level systems are valid. 

PACS numbers: 03.65.Ca, 42.50.Fx, 42.50.Ct 



I. INTRODUCTION 



The theoretical analysis of any non-trivial physical 
problem typically requires the use of approximations. A 
key approximation facilitated in most areas of physics re- 
duces the complete configuration space of the system of 
interest to a smaller set of relevant system states. In the 
theoretical description of atom-field interactions, the es- 
sential state approximation entails neglecting most of the 
bound and continuum atomic states [3, H, Q . The sem- 
inal Jaynes-Cummings-Model [1] takes this reduction to 
the extreme in that only two atomic states are retained. 
Obviously, it is essential to in detail explore the validity 
range of this reduction of the configuration space. The 
few-level approximation usually leads to theoretical pre- 
dictions that are well verified experimentally [l], Q , and 
is generally considered as understood for single-atom sys- 
tems. It fails, however, to reproduce results of quantum 
electrodynamics, where in general all possible interme- 
diate atomic states need to be considered in order to 
obtain quantitatively correct results The situation 
becomes even less clear in collective systems, where the 
individual constituents interact via the dipole-dipole in- 
teraction, despite the relevance of collectivity to many ar- 
eas of physics. Examples for such systems can be found in 
ultracold quantum gases @, trapped atoms @,[H) or solid 
state systems fol. floj. w ith applications, e.g., in quantum 
information theory [11| . 

Therefore, we discuss the validity of the few-level ap- 
proximation in dipole-dipole interacting collective sys- 
tems. For this, we study the archetype case of two dipole- 
dipole interacting atoms, see Fig. [Ha). Experiments of 
this type have become possible recently @, Q • In order to 



remain general, each atom is modelled by complete sets of 
angular momentum multiplets, as shown in Fig. [ljb). We 
find that the few-level approximation in general leads to 
incorrect predictions if it is applied to the magnetic sub- 
levels of this system. For this, we first establish a general 
statement about the system behavior under rotations of 
the atomic separation vector R. As a first conclusion 
from this result, we derive the intuitive outcome that 
the dipole-dipole induced energy shifts between collec- 
tive two-atom states are invariant under rotations of the 
separation vector R. This result can only be established 
if complete and degenerate multiplets are considered and 
dipole-dipole interactions between orthogonal transition 
dipole moments are included in the analysis. On the con- 
trary, the artificial omission of any of the Zeeman sub- 
levels of a multiplet leads to a spurious dependence of 
the energy shifts on the orientation, and thus to incor- 
rect predictions. 

For example, if in the well-known two-level approxima- 
tion only one excited state |e) and the ground state \g) 
are retained, then we recover the position-dependent en- 
ergy splitting between the entangled two-particle states 
(\e,g) ± \g,e))/V% that has previously been reported 
for a pair of two- level systems 0, E|. This geometry- 
dependence is at odds with the rotational invariance of 
the collective energy splitting expected for the degenerate 
system with all Zeeman sublevels. We thus conclude that 
the few-level approximation in general cannot be applied 
to this system. 

Our results can be generalized to more complex angu- 
lar momentum multiplets. 



II. THE MODEL 
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We describe each atom by a So «-> P\ transition shown 
in Fig. QJb) that can be found, e.g., in Ca atoms. We 
choose the z axis as the quantization axis, which is distin- 
guished by an external magnetic field that induces a Zee- 
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FIG. 1: (a) The system of interest is comprised of two 
identical atoms that are located at ri and r2, respectively. 
R = T2 — ri is the relative position of atom 2 with respect to 
atom 1. (b) Level structure of atom fi G {1, 2} which we em- 
ploy to illustrate our results. The ground state is a So singlet 
state, and the three excited levels are Zeeman sublevels of a 
Pi triplet. S is the frequency splitting of the upper levels. 



We determine the dipole moments d, — (i\d\4) via the 
Wigner-Eckart theorem [12J and find 



di = 2?e<+>, 
d 2 = Ve z , 



d 3 = -2?e<->, 



(5a) 
(5b) 



where T> is the reduced dipole matrix element and the 
circular polarization vectors are e^' = [e x ± ie y )/^/2. 

We now adapt the standard derivation of a master 
equation p| 0, Q to our multilevel system. For this, we 
assume that the radiation field is initially in the vacuum 
state denoted by g-p and suppose that the total density 
operator factorizes into a product of gp and the atomic 
density operator g at t = 0. The master equation for the 
reduced atomic density operator in Born approximation 
then takes the form 



man splitting S of the excited states. The orientation of 
R is defined relative to this quantization axis. We begin 
with the introduction of the master equation which gov- 
erns the atomic evolution of the system shown in Fig. [1] 
The internal state ji^) of atom /i is an eigenstate of 
where is the angular momentum operator of atom \i 
(/i G {1,2}). In particular, the P\ multiplet with J = 1 
corresponds to the excited states |1 M ), |2 M ) and |3 M ) with 
magnetic quantum numbers m = —1, and 1, respec- 
tively, and the So state is the ground state |4^) with 
J = m = 0. The raising and lowering operators on the 



transition of atom fj, are (i G {1, 2, 3}) 



lv)(4p 



and S. 



O) 



(1) 



The total system Hamiltonian for the two atoms and the 
radiation field is H = Ha + Hp + V , where 

3 2 

Ha = h HY. ^ S i+ S i- > Hf = Y1 fc"kaL°*« . 



i—l fi—1 

V = - d {1) ■ E{n) - d (2) ■ E(r 2 ) 



ks 



(2) 

In these equations, Ha describes the free evolution of the 
two identical atoms, Huji is the energy of state \i^) and 
we choose fkv^ = 0. Hp is the Hamiltonian of the vacuum 
field and V describes the interaction of the atoms with 
the vacuum modes in dipole approximation. The electric 
field operator E is defined as 



E(r) 



^ ik . r 

2e a v 



H.c. 



(3) 



where afc s (a ks ) are the annihilation (creation) opera- 
tors that correspond to a field mode with wave vector 
fc, polarization e ks and frequency uju, and v denotes the 
quantization volume. The electric-dipole moment oper- 
ator of atom fi is a vector operator with respect to the 
angular momentum operator of atom fi and reads 



(4) 



d t g 



■[H a ,q] 



(6) 



1 



drTr 



P^[V,U(T)[V,g F g(t-T)]U^(T)]), 



where U(r) = exp[—i(H A + Hp)r/H] and Trp() denotes 
the trace over the vacuum modes. We evaluate the inte- 
gral in Eq. ([6]) in Markov-approximation [l| and ignore all 
terms associated with the Lamb shift of the atomic levels. 
In addition, we employ the rotating-wave approximation 
and neglect anti-resonant terms that are proportional to 
SftsW and sWsW. We finally obtain 



1 2 

9tQ = -t[-Ha,0] - -AH^q] + C^g. (7) 

In this equation, the Hamiltonian Hci describes the co- 
herent part of the dipole-dipole interaction and reads 



He 



O q( 2 )cW 
s L n °i + °j - 



H.c.} 



The coefficients Qy are defined as [H, [l4| 



dj xJR) d) 



(8) 



(9) 



and the tensor \ rc 1S the real part of the tensor \ whose 
components \m f° r ^ G {!> 2, 3} are given by 



ft o 

47T£n 



9iW)hi - giiji) 



RkRi 



(10) 



Here the vector R denotes the relative coordinates of 
atom 2 with respect to atom 1 [see Fig. Q^a)], rj = koR 
and gi — (ri^ 1 +ii]~ 2 — ?/~ 3 ), 52 = (?7 _1 +3i?7~ 2 — 3?7~ 3 ) . In 
the derivation of Eq. (fTU|) , the three transition frequencies 
wi, ^2 and CJ3 have been approximated by their mean 
value lu — ck (c: speed of light) [HI]. This is justified 
since the Zeeman splitting S is much smaller than the 
optical transition frequencies a;,-. 
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The last term in Eq. {7J accounts for spontaneous emis- 
sion and reads 



1 1 

2 3 



(11) 



The total decay rate of the excited state \i) of each of the 
atoms is given by 2ji, where ji — |d.;| 2 a;Q/(67reo^c 3 ) = 
7 and we again employed the approximation LUi » loq. 
The collective decay rates result from the vacuum- 
mediated dipole-dipole coupling between the two atoms 
and are determined by 



atom 1 




R 



atom 2 



<r + -polarized 
laser field 



t 



FIG. 2: (Color online) Setup considered in Section [TTTJ which 
provides an illustration of the physical mechanisms responsi- 
ble for the breakdown of the few-level approximation. In this 
example, an external laser field is used for the sake of illustra- 
tion. Our main results starting from Section IIVI do not rely 
on external driving fields. 



r = - 

v ~ h 



(12) 



where x im = Im% is the imaginary part of the ten- 
sor x- Note that the cross terms (i ^ j) in Eqs. ([9|) 
and (|12|) represent couplings between transitions with or- 
thogonal dipole moments. If the master equation is 
transformed into the interaction picture with respect to 
Ha, terms proportional to these cross terms rotate at fre- 
quencies ±<5 or ±2(5. It follows that the parameters rib- 
and Tij are negligible if the level splitting 8 is large, i.e. 

|<5|»|ry,|r y -| (i^j). 

Next we provide explicit expressions for the coupling 
constants fiy and the decay rates in Eqs. © and 
(p~2|) . respectively. For this, it is convenient to express 
the relative position of the two atoms in spherical coor- 
dinates, 



R = R (sin cos 0, sin sin 0, cos0) . (13) 



Together with Eqs. (|10[) and ([5]) we obtain 



fi 31 = 7 _ [{rj 2 - 3) cost? - 3?/ sin 77] sin 2 0e~ 2 ^ , 

fin = 3^ [{3r/ 2 - 1 + {rj 2 - 3) cos 20) cos 77 

— 77 (1 + 3 cos 20) sin 77] , 

fi 21 = -V2cot6» fi 3 ie 10 , 

fi 22 = fin - (2 cot 2 - l)fi 3 ie 210 , 

fi 32 = -fi 2 i , fi 33 = fin , (14) 



and the collective decay rates evaluate to 



T 31 = 7— [{rj 2 - 3) sin r] + 3ricosr/] sin 2 0e~ 2 ^ . 



■_X. 



[(3t7 2 - 1 + (r/ 2 - 3) cos 20) sin 77 



+77 (1 + 3 cos 20) cos 77] , 

r 2 i = -\/2cot0r 3 ie^, 

r 22 = rn-(2cot 2 0-i)r 3 ie 2 ^, 

r 32 = -r 2i , r 33 = Tn . (15) 

A numerical study of these coupling terms can be found 
in 0. 

III. PHYSICAL MOTIVATION 



In the following Section HVi we will provide a rigorous 
treatment of the behavior of our model system under ro- 
tations of the atomic separation vector in order to study 
the geometrical properties of the different coupling terms 
in the master equation ([7]) . In order to motivate this anal- 
ysis, in this Section fllli we will discuss a simple example 
for our results. This example employs an external laser 
field driving the atoms, which is used, however, only for 
the sake of illustration. Our main results starting from 
Section HVl will not rely on external driving fields. 

To this end, we consider the geometrical setup shown 
in Fig. [5] The atoms with internal structure as in Fig. [1] 
are aligned along the y axis, and in addition to the 
model considered so far, a a + polarized laser beam with 
Rabi frequency fi^ and frequency wi, propagates in z- 
direction. In rotating-wave approximation, the atom- 
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FIG. 3: (Color online) Population in the subspace S ob- 
tained by applying the few-level approximation to the setup 
in Fig. [2] The common parameters are 6 — n/2, (f> = n/2, 
<5 = 0. Further, in (I), Ql = 2>y, R = 0.3 A , and A = 0.58 7, 
where A = ujl — wq. Curve (II) shows the case SIl = 5.4 7, 
R = 0.1A , and A = 5. 2 7. 



laser interaction Hamiltonian reads 
2 



The transition operators are defined in Eq. (JTJ) . 

Since the laser polarization is cr + , it couples only to the 
transition |3) <-> |4) in each atom. To describe this setup, 
one might be tempted to employ the usual few-level ap- 
proximation, and thus neglect the excited states |1) and 
1 2) in each atom, since they are not populated by the 
laser field. If this were correct, the seemingly relevant 
subsystem would be 

S = Span(|4,4), |3,3), |3,4), |4,3)). 

However, it is easy to prove that the state space of the 
two atoms can not be reduced to the subspace S, i.e., that 
the few-level approximation cannot be applied in its usual 
form. In order to show this, we include the atom-laser 
interaction into the master equation ([7]) and transform 
the resulting master equation in a frame rotating with 
the laser frequency. This equation is solved numerically 
with the initial condition g(t = 0) = |4,4)(4,4|, i.e. it 
is assumed that both atoms are initially in their ground 
states. 

Figure [3] shows the total population confined to the 
subspace S, 



(P S ) =Tr g(t)P s 



(16) 



where Ps is the projector onto the subspace S. It can 
easily be seen that for both sets of parameters, popula- 
tion is lost from the subspace S. Since all states but the 
excited states |1) and \2) are contained in S, it is clear 



FIG. 4: (Color online) Population of the subspace V , which 
contains all population which was lost from subspace S in 
Fig. [3] such that the population in S + V remains unity for 
all times. The parameters are as in Fig. [3] 



that it is not sufficient to take only the excited state |3) 
into account in the usual few-level approximation. 

The explanation of this outcome is straightforward. 
According to Eq. ([8]), the dipole transition |3) <-> |4) of 
one atom is coupled by the cross-coupling term ^31 to 
the |1) <-> 1 4) transition of the other atom. This cou- 
pling results in a population of state |1), even though 
the transition dipoles of the two considered transitions 
are orthogonal. Consequently, the dipole-dipole interac- 
tion between transitions with orthogonal dipole moments 
will result in the (partial) population of the states |1, 1), 
|1,3), |3, 1), 1 1 , 4) , |4, 1), although none of these states is 
directly coupled to the laser field. 

The numerical verification of these statements is shown 
in Figure 21 which depicts the population of the subspace 

V = Span(|l,l), |1,3), |3,1), |1,4), |4,1)). 

Py is the projector onto the subspace V, and the param- 
eters are the same as above. Note that we have verified 
that all population is contained in the subspace S + V , 
i.e. ( P s ) + ( Pv ) = 1 at all times. 

It is important to note that the sufficient subspace 
S + V still does not contain all possible states of the 
two atoms, because the excited state \2) of each atom is 
neglected. The justification for this is that in the chosen 
geometry, the cross-coupling terms O21, T 2 i and £1 32 , r 32 
vanish such that the transition |2) «-> |4) of one atom is 
not coupled to the transitions |1) «-> |4) and |3) <-> |4) of 
the other atom, see Eqs. (|14|) and (fT5|) . This is important 
since it demonstrates that it is also not correct to simply 
state that all atomic states have to be taken into account 
for all parameter configurations. 

The above example clearly demonstrates that the few- 
level approximation is rendered impossible by the cou- 
pling terms between transitions with orthogonal dipole 
moments. Therefore, it is the nature of the dipole-dipole 
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coupling itself which enforces that generally all Zeeman 
sublevels have to be taken into account, and not the 
polarization of the external laser fields, as one may be 
tempted to assume in the usual few-level approximation. 

A physical interpretation for the origin of the vacuum- 
induced coupling of transitions with orthogonal dipole 
moments has been given in [14]. In essence, these cou- 
plings occur if the polarization of a (virtual) photon emit- 
ted on one of the transitions in the first atom has non- 
zero projection on different dipole moments of the second 
atom. Pictorially, then the second atom cannot measure 
the polarization of the photon, and thus has finite proba- 
bility to absorb it also on transitions with dipole moments 
orthogonal to the dipole of the emitting transition. 

IV. BREAKDOWN OF THE FEW-LEVEL 
APPROXIMATION 

In this section, we return to our original setup in Fig. [TJ 
and thus drop the external driving fields employed in 
Sec. IIIII We first derive a general statement about the 
behavior of the master equation (|7|) under rotations of 
the separation vector R. On first sight, we will prove an 
obvious result: In the absence of external fields but the 
isotropic vacuum, there is no distinguished direction in 
space. Thus one expects the eigenenergies of the system 
to be invariant under rotations of R, and this is indeed 
what we find. But despite its intuitiveness, this state- 
ment needs proof, and the discussion of the proof and 
its assumptions will provide the theoretical foundation 
for our central results and physical interpretations in the 
following sections. 

A. Central theorem 

In addition to a given relative position R of the two 
atoms, we consider a different geometrical setup where 
the separation vector P is obtained from R by a rotation, 
P = lZ u (a)R. Here, TZ u (a) is an orthogonal 3x3 matrix 
that describes a rotation in the three-dimensional real 
vector space R 3 around the axis u by an angle a. Our 
aim is to show that there exists a unitary operator W 
such that 

H n (P) = WH n (R)W\ (17a) 
C y (P) Q = W [£^(R)W^gW] W\ (17b) 

where W = W u (a) is given by 

W u (a) = exp[-ia J (1) • u/h] exp[-ia J (2) ■ u/h] . (18) 

Here the operator exp[— ia J^' ■u/h] describes a rotation 
around the axis u by an angle a in the state space of 
atom [i. The notation Hq(R) and C 1 (R) means that the 
coupling constants and collective decay rates in Eqs. (O 
and (jlip have to be evaluated at R. 



We proceed with the proof of Eq. (fT"7|) . In a first step, 
we introduce the auxiliary operator Ar = WVrW^ , 
where Vr is the interaction Hamiltonian for a relative 
position of the atoms given by R, and W — W u (a) is 
defined in Eq. ([TSJ) . The evaluation of Ar involves only 
the transformation of the dipole operator of each atom. 
Since the matrix elements of vector operators transform 
like classical vectors under rotations (see, e.g., Sec. 3.10. 
in [l2|), we find 

3 

W d M Wt = J2 [di S ( 4 + H.c] , (19) 

i=i 

where di = 7?.~ 1 (a)<ii. This shows that the only differ- 
ence between the auxiliary operator Ar and Vr is that 
the dipole moments of the former are determined by di 
instead of di. 

In a second step, we employ the tensor properties of 
X to find the following expression for the parameters 
fiy(P) and Tij(P) [see Eqs. © and (fT3j)], 

m^P) = [7e- 1 (a)d i ] T x re (il) [KzH<*)dj] , (20) 

HTaiP) = [7e- 1 (a)d i ] T x lm (P) [n u \a)d*\ . (21) 

This important result shows that a rotation of the dipole 
moments di by T?." 1 ^) is formally equivalent to a rota- 
tion of R by lZ u (a) in the master equation ([7]). 

From the combination of the results obtained in step 
one and two, we conclude that the exchange of Vr by 
Ar in the integral of Eq. ([6]) is equivalent to a rotation 
of the separation vector from Rio P — lZ u (a)R, 

t 

I= J^J dTT ^{[ A *> [U(r)A R U^r),g(T)]]) (22) 
o 

= - l -[H n (P),g]+C-y(P)g, (23) 

where f = t — r and 

£(f) = t/(T)[0 F ®0(f)]C/t( T ). (24) 

Note that the equality of Eqs. ([22) and ([23]) holds under 
the same assumptions that led from Eqs. (j6]) to (JTJ). 

In the second part of the proof we evaluate the integral 
in Eq. (|2"2"|) in a different way. In the discussion follow- 
ing Eq. (dU]), we justified that £ 7 and Hq depend only 
on the mean transition frequency loq. Here we employ 
exactly the same approximation [la ] and replace the fre- 
quencies u>i appearing in U(t)ArU^{t) by ojq. Since Ha 
commutes with J^' if all frequencies Wj are replaced by 
the mean transition frequency lu q , we have [W, U] = 
and hence 

U{t)A r U\t) = WU{t)V r U\t)W^ . (25) 

It follows that the argument of the trace in Eq. ([22[) can 
be written as 

W[V R , [U{t)VrU\t),W^q{t)W]]W^ . (26) 
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In contrast to Eq. (f22|) . the double commutator contains 
now the original interaction Hamiltonian Vr that corre- 
sponds to a setting with separation vector R. We thus 
obtain 

1 = -i\WHa(,R)Wi, g] + W [C 1 {R)W ] qW) . (27) 

Finally, the comparison of Eqs. (|27|) and (|23|) establishes 
Eq. (fT7|) which concludes the proof. 

Note that throughout this proof, we have not made 
reference to the specific type of the Zeeman sublevels 
employed in our example shown in Fig. [TJ Therefore, 
the central theorem holds for transitions between states 
with arbitrary angular momentum structure, as long as 
complete multiplcts are considered. 

B. Diagonalization of Hn 

We now turn to the discussion of Eq. (flT]) . which will 
lead to our central results. The Hamiltonian Hq de- 
scribes the coherent part of the dipole-dipole interaction 
between the atoms. From Eq. (|17a|) , it is immediately 
clear that the eigenvalues of Hq depend only on the in- 
teratomic distance, but not on the orientation of the sep- 
aration vector R. The reason is that the spectrum of two 
operators, which are related by a unitary transformation, 
is identical. In our case, the Hamiltonian Hq(R) and 
Hq (P) for different orientations R and P are related by 
the unitary transformation W, and since P is obtained 
from R by an arbitrary rotation, the eigenvalues of Hq 
are identical for any orientation. 

Next we re-obtain this result in a more explicit way 
and derive symbolic expressions for the eigenvalues and 
eigenstates of Hq. This Hamiltonian can be written as 

3 

Hq = ]T ( l fl *X fl il + k>fal ) . (28) 

ij = l 

where the symmetric and antisymmetric states are de- 
fined as 

|« i ) = (|i,4) + |4,i))/V2, (29a) 
\ ai ) = (\i,A)-\A,i))/V2, (29b) 

and \i,j) = \ii) <8> | j 2 ) ■ Since all matrix elements 
(sj|i?n|aj) of Hq between a symmetric and an antisym- 
metric state vanish, the set of eigenstates decomposes 
into a symmetric subspace S and an antisymmetric sub- 
space A. The matrix elements of Hq in the subspace iS 
spanned by the symmetric states {|si), | S2) , \ss)} are 

(fin H|i ^31 \ 
n 21 n 22 n* 32 , (30) 
^31 f232 ^33 J 

and the representation of Hq in the subspace A spanned 
by the antisymmetric states {|ai), \a 2 ), {0,3)} is given by 



[Hq] a — — [Hq] s . Note that the collective ground state 
|4,4) and the states \i,j) £ {1,2,3}) where each 
atom is in an excited state are not influenced by the 
dipole-dipole interaction and thus not part of the expan- 
sion ([28|). 

In Section HV Al we have derived a general relation be- 
tween any two orientations of the interatomic distance 
vector. In order to apply this result, we define the vector 
R z to be parallel to the z axis, i.e. R z = Re z . This cor- 
responds to the choice 9 = in Eq. (fT3|) . Any separation 
vector P can then be obtained from R z as P = lZ u (a)R z 
by a suitable choice of the rotation axis it and the angle 
a. 

We then proceed with the diagonalization of the Hamil- 
tonian Hq(R z ) with atomic separation vector R z . The 
explicit calculation of the coupling constants f2y shows 
that the off-diagonal elements in Eq. (|30|) vanish if the 
atoms are aligned along the z axis, see Eqs. (p~4|) and (fl5|) 
with 9 = 0. It follows that the Hamiltonian Hq(R z ) is al- 
ready diagonalized by the symmetric and antisymmetric 
states Eq. ([29]) . and the eigenvalues of [iJsi] 5 and [i?h] 
are given by Xf = —hQa(R z ) and = Mla(R z ), re- 
spectively. 

According to Eq. (|17a|) , the Hamiltonian Hq(P) is the 
unitary transform of Hq(R z ) by W. The normalized 
eigenstates of Hq(P) are thus determined by W|sj) and 
W|aj), and their eigenvalues are again Xf and Xf , re- 
spectively. Since the orientation of P is arbitrary, the 
eigenvalues of Hq(P) depend only on the interatomic 
distance \P\ = \R Z \ = R, but not on the orientation of 
the separation vector. 

Thus, it follows from our theorem in Sec. lIV Al that the 
eigenvalues of Hq (P) are invariant under rotation of the 
interatomic distance vector. 



C. Diagonalization of Ha + Hn 

An additional conclusion can be drawn from Eq. (|17p 
if the operator Ha commutes with the transformation 
W = W u (a), i.e., 

[H A ,W]=0. (31) 

Then, Eq. (|17ap implies that Ha + Hq(P) is the unitary 
transform of Ha + Hq(R) by W. A straightforward re- 
alization of this is the case of vanishing Zeeman splitting 
5, in which the relation holds for an arbitrary orientation 
of P. Then, the energy levels of the full system Hamil- 
tonian Ha + Hq do not depend on the orientation of the 
separation vector. 

This result can be understood as follows. In the ab- 
sence of a magnetic field (6 = 0), there is no distinguished 
direction in space. Since the vacuum is isotropic in free 
space, one expects that the energy levels of the system 
are invariant under rotations of the separation vector R. 

By contrast, the application of a magnetic field in z 
direction breaks the full rotational symmetry. For S =/= 0, 
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the atomic Hamiltonian Ha only commutes with trans- 
formations W u (a) that correspond to a rotation of the 
separation vector around the z axis, u — e z . If we ex- 
press the atomic separation vector in terms of spherical 
coordinates as in Eq. (fT3|) , this means that the eigen- 
values of the full system Hamiltonian Ha + Hn do only 
depend on the interatomic distance R and the angle 9, 
but not on the angle 4>. This result reflects the symmetry 
of our system with respect to rotations around the z axis. 



D. Unitary equivalence of time evolution in 
different orientations 



If the operator Ha commutes with the transformation 
W = W u (a), another conclusion can be drawn. Then, 
the result in Eq. (JT7J) implies that the density operator 
Wq(R)W* obeys the same master equation than g(P) 
for P = lZ u (a)R. It follows that P is the unitary trans- 
form of g(R) by W, i.e. 



B (P) = Wg(R)W* . 



(32) 



As discussed in Sec. IIV C[ the free atomic Hamiltonian 
Ha commutes with W u (a) for an arbitrary choice of the 
rotation axis u and angle a if the Zeeman splitting 5 
vanishes. 

We thus conclude that it suffices to determine the so- 
lution of the master equation ([7]) for only one particular 
geometry if 6 = 0. Any other solution can then be gener- 
ated simply by applying the transformation W = W u (a) 
with suitable values of u and a to the solution for the 
particular geometry. 



E. Establishment of the breakdown 



In Sees. IIV BIIIV Dl we presented several results con- 
cerning the energy levels and the time evolution of the 
two dipole-dipole interacting atoms that are based on the 
central theorem in Eq. (fT7|) . However, this theorem can 
only be established if each atom is modelled by complete 
sets of angular momentum multiplets, and represents the 
reference case that corresponds to results which can be 
expected in an experiment. If any of the Zeeman sub- 
levels of the Pi triplet are neglected, the unitary operator 
W does not exist since it is impossible to define an angu- 
lar momentum or vector operator in a state space where 
magnetic sub-levels have been removed artificially. In this 
case, the central statement cannot be applied. Still, the 
system can be solved without the help of the theorem. 
The breakdown of the few-level approximation for collec- 
tive systems is then established by noting that the results 
for systems with artificially reduced state space fail to re- 
cover the results derived in Sees. IIV BIIIV Dl for the full 
system. 

In order to illustrate this point in more detail, we con- 
sider the system in Fig. [1] and assume that the excited 
states of each atom are degenerate (8 — 0). According to 



our findings in Sees. lIVBl and UVCI the energy levels of 
the complete system depend on the length of the separa- 
tion vector i?, but not on its orientation. In contrast, the 
omission of any of the Zeeman sublevels leads to a spu- 
rious dependence of the energy levels on the orientation, 
and thus to incorrect predictions. 

For example, if the excited states |1) and |3) in each 
atom are omitted, the level scheme in Fig.QJb) reduces to 
an effective two-level system comprised of the states |2) 
and 1 4). The collective two- atom system is then described 
by the ground state |4, 4), the excited state |2, 2) and the 
symmetric and antisymmetric states |sa) and 1^2)- The 
frequency splitting between the states and \a,2) is 
given by 2| O22 1 , where 



^22 = ^[/iW-cos 2 ^)/^)] 



and 



Mv) 



1 1 

77 T 

1 _3 

77 if 



COS 7/ 



cos 77 



_3 

if 2 



sin 77 . 



sin 77 . 



(33) 

(34) 
(35) 



Since the second term in Eq. (|33|) is proportional to 
cos 2 9, the energy levels of the artificially created two- 
level system strongly depend on the orientation of the 
separation vector R. This is at variance with our finding 
in Sec. IIV CI where we have shown that the energy lev- 
els do not depend on the orientation of the vector R if 
each atom consists of complete and degenerate Zeeman 
multiplets. We thus conclude that all Zeeman sublevels 
generally have to be taken into account. 

Since the validity of the central theorem Eq. (fT7|) is not 
restricted to the So <-> Pi transition discussed so far, it 
follows that the breakdown of the few-level approxima- 
tion can be established for transitions between arbitrary 
angular momentum multiplets. 

The intuitive explanation of the breakdown has already 
been hinted at in Sec. Mil For a more formal discus- 
sion, we return to the matrix representation of [Hq] s in 
Eq. (|3"D|) . The diagonal elements proportional to f2,-, ac- 
count for the coherent interaction between a dipole of 
one of the atoms and the corresponding dipole of the 
other atom. By contrast, the off-diagonal terms pro- 
portional to Qij with i ^ j arise from the vacuum- 
mediated interaction between orthogonal dipoles of dif- 
ferent atoms [l3|, EH • It is the presence of these terms 
that renders the simplification of the atomic level scheme 
impossible since they couple an excited state \i) of one 
atom to a different excited state \j) (i 7^ j) of the other 
atom. A similar argument applies to the collective decay 
rates T+j appearing in C 7 g. Thus, if any Zeeman sublevel 
of the excited state multiplet is artificially removed, then 
some of these vacuum-induced couplings f2y- with i ^ j 
are neglected, which leads to incorrect results. Now, it is 
also apparent why the breakdown of the few-level approx- 
imation appears exclusively in collective systems. For 
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single atoms in free space, a coupling of orthogonal tran- 
sition dipole moments via the vacuum is impossible. 

F. Recovery of the few-level approximation in 
special geometries 

The identification of the vacuum-induced couplings f2y 
and Tij between orthogonal transition dipole moments as 
the cause of the breakdown enables one to conjecture that 
few-level approximations are justified for particular geo- 
metrical setups, where some or all of the cross-coupling 
terms vanish. 

For example, we mentioned earlier that all cross- 
coupling terms vanish if the atoms are aligned along the 
z axis. This corresponds to the case 9 = in Eqs. (fla - 
p's]) . Then, the Sq <-> Pi transition may be reduced to a 
two-level system, formed by an arbitrary sublevel of the 
Pi triplet and the ground state Sq. 

As a second example, we assume the atoms to be 
aligned in the x-y-plane, i.e., 9 = ir/2 in Eq. (| 13|) . Then 
the terms f2 21 , r 21 and f2 32 , r 32 vanish, see Eqs. (fTi)) - 
(dnj). In effect, the excited state \2) may be disregarded 
such that the atomic level scheme simplifies to a V-system 
formed by the states |1) and |3) of the Pi multiplet and 
the ground state Sq. 

Note that the cross-coupling terms also become irrele- 
vant in the special case \S\ 3> |ry| (i ^ j), see our 
discussion below Eq. (fT^l) . 

V. DISCUSSION AND SUMMARY 

Throughout this article, we have studied the properties 
of various parts of the system Hamiltonian as well as the 
full density operator under rotations of the interatomic 
distance vector. This discussion was based on a general 
theorem in Sec. IIV Al which relates the system proper- 
ties for different orientations of the interatomic distance 
vector. 

First, we have discussed the Hamiltonian Hq, which 
describes the coherent coupling between different tran- 
sitions in the two atoms induced by the vacuum field. 
Armed with our main theorem, it is possible to first di- 
agonalize Hq in a special geometry, where the eigenvec- 
tors and eigenenergies assume a particularly simple form. 
The eigenvectors and eigenenergies for an arbitrary sys- 
tem geometry are then derived via the theorem. Our 
main result of Sec. IIV Bl is that the eigenvalues of Hq 
are invariant under rotation of the interatomic distance 
vector. 

In a second step, we have studied the eigenenergies 
of the full system Hamiltonian Ha + Hq , which in gen- 
eral are not invariant under rotation of the interatomic 
distance vector. The invariance, however, is recovered 



if Ha commutes with the transformation W — W u (a), 
which is given in explicit form as a result of our theorem. 
Most importantly, this additional condition is fulfilled for 
a degenerate excited state multiplet, i.e., if the Zeeman 
splitting S vanishes. Then, there is no preferred direction 
in space, such that the invariance of the eigenenergies, 
which are observables, can be expected. 

We then conclude the breakdown of the few-level ap- 
proximation in Sec. IIV fi] since our results of the previous 
sections are violated if any of the excited state multiplet 
sublevels are artificially removed. Possible consequences 
are, for example, a spurious dependence of the eigenener- 
gies on the orientation of the interatomic distance vector, 
and thus of all observables that depend on the transition 
frequencies among the various eigenstates of the system. 
In experiments, in addition, a loss of population from 
the subspace considered in the few-level approximation 
would be observed. Our proof can be generalized to tran- 
sitions between arbitrary angular momentum multiplets. 

We have identified the vacuum-induced dipole-dipole 
coupling between transitions with orthogonal dipole mo- 
ments as the origin of the breakdown. On the one hand, 
this explains why the breakdown exclusively occurs in 
collective systems, since such orthogonal couplings are 
impossible in single atoms in free space. On the other 
hand, the interpretation enables one to identify special 
geometries where some of the Zeeman sublevels can be 
omitted. This also allows to connect our results to previ- 
ous studies involving dipole-dipole interacting few-level 
systems. In these studies involving the few-level approx- 
imation, typically a very special geometry was chosen, 
e.g., with atomic separation vector and transition dipole 
moments orthogonal or parallel to each other. These re- 
sults remain valid if a geometry can be found such that 
the full Zeeman sublevel scheme reduces to the chosen 
level scheme as discussed in Sec. IIV Fl It should be noted, 
however, that there are physical realizations of interest 
which in general do not allow for a particular system 
geometry that leads to the validity of a few-level approx- 
imation, such as quantum gases. 

Finally, on a more technical side, our results can also 
be applied to considerably simplify the computational 
effort required for the treatment of such dipole-dipole in- 
teracting multilevel systems with arbitrary alignment of 
the two atoms. First, our theorem both allows for a con- 
venient evaluation of eigenvalues and eigenenergies for 
arbitrary orientations of the interatomic distance vector 
based on the results found in a single, special alignment. 
Second, we have found in Sec. IIV Dl that for the degener- 
ate system, the density matrices for different orientations 
are related to each other by the unitary transformation 
W defined in our theorem. Thus the solution for any ori- 
entation can be obtained from a single time integration 
simply by applying this transformation. 
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